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From neutron diftraction measurements on a quasi-lD Ising-like Co^+ spin compound 
BaCo2V2 08, we observed an appearance of a novel type of incommensurate ordering in magnetic 
fields. This ordering is essentially different from the Neel-type ordering, which is expected for 
the classical system, and is caused by quantum fluctuation inherent in the quantum spin chain. 
A Tomonaga-Luttinger liquid (TLL) nature characteristic of the gapless quantum ID system is 
responsible for the realization of the incommensurate ordering. 
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The quasi-one-dimensional (ID) antiferromagnets, in 
which spin chains are coupled only weakly by interchain 
interactions, involve significant quantum fluctuation and 
often display exotic behavior. As expected for the classi- 
cal models, the antiferromagnet in three-dimension typi- 
cally develops a usual Neel order with antiparallel align- 
ments of neighboring spins when it is cooled. Owing to 
the strong quantum fluctuation, however, qualitatively 
different situations ft'om the classical models can appear 
for the quasi-lD antiferromagnets. Several exotic states 
have been discovered in the quasi- ID antiferromagnets, 
such as the valence-bond-solid ground state |lj or the 
spin-Peierls state [3]. In this letter, we show that in the 
system with Ising-like anisotropy in magnetic fields, the 
symmetry breaking for the long range order takes place 
in a very unusual way. We have found a novel type of 
density- wave-like incommensurate ordering in a quasi- ID 
Ising-like Co^"'" spin compound BaCo2V2 08. This curi- 
ous ordering originates from instability of the peculiar 
quantum critical nature, characterized by an incommen- 
surate modulation wave vector in one-dimensional chain, 
that has no classical analog. 

The aniferromagnetic chain in magnetic fields applied 
along the longitudinal z-direction is described by the fol- 
lowing XX Z Hamiltonian: 

i 

-gfj'B^Si.zH (1) 

i 

where, J (>0) is an antiferromagnetic exchange constant, 
e an anisotropic parameter, g a value, the Bohr mag- 
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FIG. 1: Spin correlations in the TLL phase of the S = 
1/2 ID XX Z antiferromagnet. (a) The staggered correla- 
tion of the transverse spin components, which is expressed 
as < SqS^ > ~(— 1)'V~''''. (b) The incommensurate cor- 
relation of the longitudinal spin component, expressed as 
< SqS^ > — m^~cos(2fcFr)r~''^ . r is a distance between the 
spins. T^x and r/z are Tomonaga-Luttinger exponents, which 
satisfy a relation ri^riz ~ 1. 



neton and H the magnetic field. The system with e < 1, 
e = 1 and e > 1 corresponds to the Ising-like, Heisenberg 
and Xy-like antiferromagnetic chain, respectively. The 
quantum fluctuation, stemming from non-commutative 
properties of spin operators, plays a crucial role in de- 
termining its ground state properties. Indeed, the ex- 
act quantum mechanical ground state of the Heisenberg 
chain with spin S = 1/2, found by Bethe in 1931 3*1, is 
a spin liquid state with no order, showing that the quan- 
tum fluctuation in this case disturbs the system in taking 
a long range ordered state even at absolute zero Kelvin. 
After a theoretical flnding that the spin chain with S — 
1/2 can be represented by a pseudo-fermion modelfl, [H], 
an important conclusion has been extracted ^, '7^7 i-e. 
the spin liquid state belongs to a universality class called 
a Tomonaga-Luttinger liquid (TLL). The TLL includes 
most ID quantum systems having a gapless elementary 
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FIG. 2: Phase diagram of BaCo2V208 in magnetic fields ap- 
plied along the c-axis. Neutron scattering measurements are 
conducted in the region, shown by allows. 



excitation with a linear energy dispersion, such as a lin- 
ear chain conducting electron system. In general, a TLL 
has no long range order even at zero Kelvin, but it is 
in a quantum critical state with intrachain correlations 
of algebraic decay It\. A distinctive feature of the TLL 
in the spin chain is that two kinds of the correlation de- 
velop there as shown in Fig. 1 [1, 0, B @- One is the 
staggered correlation of the transverse spin component 
perpendicular to the chain direction. The other is the in- 
commensurate correlation of the longitudinal component 
parallel to the chain. One may associate the former with 
the usual Neel order, but the later has no classical analog 
and is peculiar to the quantum spin chain. 

The TLL nature which appears in the S = 1/2 Ising- 
like antiferromagnetic spin chain, brings the curious or- 
dering found in this study. The Ising-like chain shows 
an interesting quantum phase transition in magnetic 
fields @, • Different from the Hcisenberg system, the 
ground state of the Ising-like chain has a long range order 
at zero magnetic field, because the Ising-like anisotropy 
stabilizes the Neel state [ll|. The external field, however, 
induces the strong quantum fluctuation, that drives the 
system into the TLL phase above a certain critical field. 
The idea underlying our results concerns with the real 
quasi-lD compounds, in which interactions between the 
chains inevitably exist [l^, [l^ [ij, [3 ■ The interchain in- 
teractions make the TLL phase unstable, and can thereby 
lead the system into the long range ordering at a finite 
temperature. On this situation, between two kinds of the 
correlation mentioned above, the most dominant one will 
grow rapidly and build up the order. An important point 
for the Ising-like chain is that the incommensurate cor- 
relation is enhanced above the critical field, in contrast 
to the Heisenberg system, for which the staggered corre- 
lation is always dominant in a magnetic field p^ . Thus, 
we propose that a density-wave-like incommensurate or- 
dering appears for the quasi ID Ising-like antiferromag- 
net in the field-induced region, when the temperature is 



lowered enotigh to make the interchain interaction rele- 
vant [H, 0, In this ordered state, the spins ahgn to 
be collinear along the chain direction with modulation of 
those amplitude, characterized by the incommensurate 
wave number 2kp. This ordering is essentially different 
from that expected for the classical Ising-like antiferro- 
magnet with no quantum fluctuation. The field-induced 
transition in the classical case is a spin-flop type, and 
the Neel order of the transverse component, which re- 
lates to the breaking of rotational symmetry around the 
field, appears in the field-induced region as well as the 
Heisenberg antiferromagnct [17]. On the other hand, the 
incommensurate ordering discussed here occurs as a con- 
sequence of the breaking of quasi-continuous translation 
symmetry. According to the pseudo-fermion model, this 
incommensurate ordering is regarded as a charge den- 
sity wave (CDW) ordering for the pseudo-fermion. The 
modulation wave number 2kp of this incommensurate or- 
dering can be easily tuned by applying magnetic fields, 
since it is determined by the magnetization per spin m 
as 2kF = 7r(l - 2m) This means that the periodic- 
ity of the incommensurate structure varies continuously 
with the field strength. This property is also different 
from that of the density-wave ordering in the conducting 
electron systems, of which 2kF depends on the electron 
occupation number and is therefore little affected by ex- 
ternal perturbations. The density-wave ordering in the 
conducting electron systems is arisen by the nesting of 
Fermi surface, whereas the incommensurate ordering in 
the quasi ID Ising-like XXZ case is entirely due to the 
quantum fluctuation inherent in the system. 

In order to find this curious ordering, neutron diffrac- 
tion measurements, which provide direct information 
about the spin structure, are particularly useful. For 
this measurement, we adopt the Co^"*" spin system 
BaCo2V208, in which magnetic Co-0 chains are running 
along the crystallographic c-axis [11] . Recently, this com- 
pound was revealed to be a good realization of quasi-lD 
5 = 1/2 Ising-like antiferromagnct with the transition 
field He ~ 3.9 T [19, 20] that can be achieved in current 
neutron facilities. The field-temperature phase diagram 
of BaCo2V208, obtained from our thermodynamic mea- 
surements [H, [l^, are depicted in Fig. 2. BaCo2V208 
undergoes Neel ordering at Tn = 5.4 K at zero magnetic 
field, but the ordered temperature is rapidly lowered by 
the external field along the chain, which corresponds the 
easy axis direction [l^ . The suppression of the Nccl order 
by magnetic fields can be understood by an appearance of 
the TLL nature in the field, which was mentioned for the 
Ising-like chain before (20| . However, we recently found 
that at very low temperatures below 1.8 K, another or- 
dered phase emerges in the field-induced region above He 
~ 3.9 T [l^. In the ordered phase in the field- induced 
region, we expect a realization of the incommensurate 
spin structure. 

We have performed the neutron diffraction measure- 
ments in the following condition. A single crystal of 
BaCo2V208, having a shape of a plate with 6 x 6 x 26 
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FIG. 3: Magnetic field dependence of neutron difi'raction pro- 
files of (4, 0, I) scan measured at temperature T — 0.85 K 
(a), 1.5 K (b) and 2 K (c). Solid lines are the results of fits 
to Gaussian functions. 



mm'^ and a weight of 1.5 g, was used for neutron elastic 
scattering measurements. A space group of BaCo2V2 08 
is 141/acd. The neutron scattering measurements were 
carried out with the thermal neutron triple-axis spec- 
trometer TAS-2 at the JRR-3 reactor at Japan Atomic 
Energy Agency. Incident neutrons are monochromatized 
by the (0 2) reflection of pyrolytic graphite (PG) crys- 
tal, and contamination from higher-order beams was ef- 
fectively eliminated using a PG filter. The measurements 
were conducted in horizontal fields up to 5 T. A split- 
pair superconducting magnet manufactured by Oxford 
Instruments, UK, was used for the field generation. The 
fixed incident neutron energy was 13.7 meV. The hori- 
zontal collimator sequence was guide-80'-80'-open. The 
cooling of the sample was achieved by a '^He-'^IIe dilution 
refrigerator. 

Now let us turn to our neutron diffraction investigation 
of the spin structure in BaCo2V2 08- The measurement 
was conducted in magnetic fields up to 5 T along the c- 
axis and at temperature below 2 K. Figure 3 shows the 
field dependence of the scan profile of the (4 /) reflection 
in the fields. In neutron scattering measurements, only 
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FIG. 4: Magnetic field dependence of the peak position of 
the observed neutron scan profile. The data are extracted 
from the least-square fits to the neutron scan profile. Inset 
shows the field dependence of a normalized incommensurate 
modulation 2kF/n. The theoretical curve is obtained by the 
calculation based on the Bethe anzats exact theory. 



the component of the magnetization that is perpendic- 
ular to the scattering vector contributes to the scatter- 
ing intensity ^27\. Thus, to detect the magnetic ordered 
component parallel to the chain direction, we adopt the 
Z-scan around (4 3). In fact, other reflections such as (4 
1) are more suitable to detect the magnetic component 
along the c axis. However, these reflections do not satisfy 
the configurational condition of our measurement system, 
which is restricted by the windows of the sprit-pair mag- 
net. At zero magnetic field, a peak at (4 3), which 
corresponds to the Neel order of the magnetic moments 
along the chain, is observed. According to the extinction 
rule, the (4 3) nuclear reflection in 141/acd space group 
is prohibited. Thus, the observed (4 3) peak is purely 
magnetic, and we confirmed a disappearance of the peak 
above Tn = 5.4 K. The (4 3) peak gradually diminishes 
with increasing the field up to 3.75 T, and then a sudden 
change of the scan profile occurs around the transition 
field He- In the field- induced region above He, two peaks 
at positions incommensurate with the underlying lattice 
appear at temperatures below 1.5 K. The sudden change 
of the scan profile reflects the fact that the transition at 
He is weakly first order as sug gest ed from our previous 
thermodynamic measurements [15|. A tiny peak at (4 
3) remains in the field-induced region, but its origin is not 
clear at the moment. The positions of the peaks are plot- 
ted in Fig. 3. The difference between the peak positions 
in field ascending and descending processes is small. The 
incommensurate peaks shift continuously in such a way 
that a distance between the two peaks increases with in- 
creasing the field. Slight temperature dependence of the 
peak positions is found for H > He- The line width for 
all the peaks, observed in the ordered region, is within 
a resolution limit of our apparatus, and we confirm that 
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the incommensurate peak, observed from the /i-scan, is 
also resolution-limited. The resolution-limited peaks cor- 
respond to the development of the long range order. The 
scan profiles, observed at 4.5 T, show that the Bragg 
peaks change to broad defuse ones at 1.5 K, that is just 
above the ordering temperature, and then disappear with 
increasing the temperature. Our observation unambigu- 
ously demonstrates a realization of the incommensurate 
ordering in the field-induced phase of BaCo2V208. In our 
experimental accuracy, no peak is found in the Z-scan pro- 
file around (0 1) at which the neutron cross section in- 
cludes solely transverse component of magnetic moment. 
The result indicates the density-wave-like ordering with 
collinear spin alignments along the chain directions in the 
field-induced phase. Satellite reflections, coming from a 
higher harmonic Fourier components of the incommen- 
surate modulation for the ordered structure, are not ob- 
served in the scan profile of the (4 /) refiection. Thus, 
the incommensurate modulation is suggested not to be 
square-wave like but is close to proper sinusoidal. In the 
inset of Fig. 4, we plot the field dependence of a normal- 
ized incommensurate modulation 2k-p /t^ along the c-axis, 
which is given by a relation ^k^jix ~ \ - 8. Taking into 
account the fact that four Co^+ ions are included in a 
chemical unit of BaCo2V2 08 along the chain, the 5 is 
obtained from 1/8 of a distance between two incommen- 



surate peaks. As anticipated, the Ik^j-K continuously 
decreases with increasing the field above H^. The ex- 
perimental result slightly deviates from the theoretical 
prediction Ik-p = 7r(l- 2to) Q with increasing field as 
shown in the inset of Fig. 4. The theoretical curve is 
obtained by the calculation based on the Bethe ansatz 
exact theory with the parameters J/Zcb = 65 K, e = 0.46 
and g = 6.2, which are estimated from the magnetization 
curve [20[ . The observed two incommensurate peaks re- 
fiect an existence of two kinds of domain for the ordered 
structure. The left and right peaks correspond to the 
domain with the modulation wave vector pointing par- 
allel and antiparalell to the field direction in the chain, 
respectively. The difference of the intensity between the 
two peaks probably originates from the imbalance in two 
domains caused by the external field. 

The authors are grateful to Dr. T. Suzuki for valu- 
able discussions. This work was partly supported by 
Grants-in-Aid for Young Scientists (B) (No. 18740183, 
No. 1870230 and No. 19740215), for Scientific Research 
on Priority Areas "High Field Spin Science in 100 T" 
(No. 451) and "Novel State of Matter Induced by Frus- 
tration" (No. 19052004 and No. 19052008) and for Sci- 
entific Research (B) (No. 20340089) from the Ministry 
of Education, Science, Sports, Culture and Technology 
(MEXT) of Japan. 



[1] I. Affleck, T. Kennedy, E.H. Lieb, and H. Tasaki, Phys. 

Rev. Lett 59 799 (1987). 
[2] E. Pytte, Phys. Rev. B 10 4637 (1974). 
[3] H. Bethe, Z. Physik. 71 205 (1931). 

[4] E. Leib, T. Shultz and D. Mattis, Ann. Phys 16 407 
(1961). 

[5] L. Bulaevskii, Soviet Phys.-JETP 16 685 (1963). 
[6] F.D.M. Haldane, Phys. Rev. Lett. 45 1358 (1980). 
[7] A. Luther and I. Paschel, Phys. Rev B 12 3908 (1975). 
[8] N. Ishimura and H. Shiba, Prog. Theor. Phys. 64 479 
(1980). 

[9] G. Muller, H. Thomas, H. Beck and J.C. Bonner, Phys. 
Rev. B 24 1429 (1981). 
[10] C.N. Yang and CP. Yang, Phys. Rev. B 150 321, 327 
(1966). 

[11] R. Orbach, Phys. Rev. 112 309 (1958). 
[12] T. Sakai, Phys. Rev. B 58 6268 (1998). 
[13] T. Suzuki, N. Kawashima and K. Okunishi, J. Phys. Soc. 

Jpn. 76 123707 (2007). 
[14] K. Okunishi and T. Suzuki, Phys. Rev. B 76 224411 



(2007). 

[15] S. Kimura, T. Takeuchi, K. Okunishi, M. Hagiwara, Z. 

He, K. Kindo, T. Taniyama and M. Itoh, Phys. Rev. Lett. 

100 057202 (2008). 
[16] N. M. Bogoliubov, A.G. Izergin and V.E. Korepin, Nucl. 

Phys. B 275 687 (1986). 
[17] M.E. Fisher and D.R. Nelson, Phys. Rev. Lett. 32 1350 

(1974). 

[18] R. Wichmann and Hk. Muller-Buschbaum, Z. Anorg. 

Allg. Chem. 534 153 (1986). 
[19] Z. He, T. Taniyama, T. Kyomen and M. Itoh, Phys. Rev 

B 72 172403 (2005). 
[20] S. Kimura, H. Yashiro, K. Okunishi, M. Hagiwara, Z. He, 

K. Kindo, T. Taniyama and M. Itoh, Phys. Rev. Lett. 99 

087602 (2007). 

[21] G. L. Squires, Introduction to the theory of thermal 
neutron scattering (Dover publications inc.. New York, 
1978). 



